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£2 ' Abstract 

o 

ly-^ , In this article we seek stability estimates in the inverse problem of de- 

termining the potential or the velocity in a wave equation in an anisotropic 

^^ , medium from measured Neumann boundary observations. This information 

is enclosed in the dynamical Dirichlet-to-Neumann map associated to the 
wave equation. We prove in dimension n > 2 that the knowledge of the 
Dirichlet-to-Neumann map for the wave equation uniquely determines the 
electric potential and we prove Holder-type stability in determining the po- 
tential. We prove similar results for the determination of velocities close 
to 1. 

Keywords: Stability estimates, Hyperbolic inverse problem, Dirichlet-to- 
Neumann map. 
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1 Introduction 

In this paper, we are interested in the following inverse boundary value problem: 
on a Riemannian manifold with boundary, determine the potential or the velocity 
— i.e. the conformal factor within a conformal class of metrics — in a wave 
equation from the vibrations measured at the boundary. Let (M, g) be a compact 
Riemannian manifold with boundary dA4. All manifolds will be assumed smooth 
(which means C°°) and oriented. We denote by A g the Laplace-Beltrami operator 
associated to the metric g. In local coordinates, 

n 

g( x ) = ^2 gjk(x)dxj (8) dx fc , 
j,k=l 



A g is given by 

A P 



Vdetg^dxj V " gS dx k ) ' 



Here {g^ k ) is the inverse of the metric g and det g = det(gjfc). Let us consider the 
following initial boundary value problem for the wave equation with bounded (real 



valued) electric potential q G L°° (M ) 

' (a t 2 - A g + q(x)) u = 0, in(Q,T)xM, 

u(0,-) = 0, d t u(0,-) = in M, (1.2) 

u = /, on(0,T)xd.A/i, 

where / € -ff 1 ((0, T) x dM). Denote by v = u(x) the outer normal to dM at 

n 

x G dM, normalized so that Y^ Z^ ky j v k = 1- We may define the dynamical 

i,fc=i 
Dirichlet-to-Neumann map A gi q by 



O'U 



W= £■*«** 



j,fc=l 



(1.3) 

(o,r)xax 



It is clear that one cannot hope to uniquely determine the metric g = (gjk) from 
the knowledge of the Dirichlet-to-Neumann map A g(? . As was noted in |[39l . the 
Dirichlet-to-Neumann map is invariant under a gauge transformation of the met- 
ric g. Namely, given a diffeomorphism \& : M — > M. such that ^lax = Id one 
has A^r* gi q = A gi q where ^*g denotes the pullback of the metric g under fy. The 
inverse problem should therefore be formulated modulo the natural gauge invari- 
ance. Nevertheless, when the problem is restricted to a conformal class of metrics, 
there is no such gauge invariance and the inverse problem now takes the form: 
knowing A cg g , can one determine the conformal factor c and the potential q? 

Belishev and Kurylev gave an affirmative answer in [5] to the general problem 
of finding a smooth metric from the Dirichlet-to-Neumann map. Their approach 
is based on the boundary control method introduced by Belishev [4] and uses in 
an essential way an unique continuation property. Unfortunately it seems unlikely 
that this method would provide stability estimates even under geometric and topo- 
logical restrictions. Their method also solves the problem of recovering g through 
boundary spectral data. The boundary control method gave rise to several refine- 
ments of the results of (H: one can cite for instance |[30ll . |29l and HI. 

In this paper, the inverse problem under consideration is whether the knowl- 
edge of the Dirichlet-to-Neumann map A gj q on the boundary uniquely determines 
the electric potential q (with a fixed metric g) and whether the knowledge of the 
Dirichlet-to-Neumann map A g = A g o uniquely determines the conformal factor 
of the metric g within a conformal class. From the physical viewpoint, our inverse 
problem consists in determining the properties (e.g. a dispersion term) of an in- 
homogeneous medium by probing it with disturbances generated on the boundary. 



The data are responses of the medium to these disturbances which are measured 
on the boundary, and the goal is to recover the potential q(x) and the velocity c(x) 
which describes the property of the medium. Here we assume that the medium is 
quiet initially, and / is a disturbance which is used to probe the medium. Roughly 
speaking, the data is d v u measured on the boundary for different choices of /. 



In the Euclidian case (g = e) Rakesh and Symes [35], [34] used complex ge- 
ometrical optics solutions concentrating near lines with any direction oj e S n_1 to 
prove that A e>q determines q(x) uniquely in the wave equation. In IJ351 . h. e ^ q gives 
equivalent information to the responses on the whole boundary for all the possible 
input disturbances. Ramm and Sjostrand [36] extended the results in [35] to the 
case of a potential q depending both on space x and time t. Isakov |[25l considered 
the simultaneous determination of a potential and a damping coefficient. A key 
ingredient in the existing results, is the construction of complex geometric optics 
solutions of the wave equation in the Euclidian case, concentrated along a line, and 
the relationship between the hyperbolic Dirichlet-to-Neumann map and the X-ray 
transform plays a crucial role. 



Regarding stability estimates, Sun B2II established in the Euclidean case sta- 
bility estimates for potentials from the Dirichlet-to-Neumann map. In 091 and 
fill Stefanov and Uhlmann considered the inverse problem of determining a Rie- 
mannian metric on a Riemannian manifold with boundary from the hyperbolic 
Dirichlet-to-Neumann map associated to solutions of the wave equation {df — 



A g )u = 0. A Holder type of conditional stability estimate was proven in Il39l 
for metrics close enough to the Euclidean metric in C k , k > 1 or for generic simple 
metrics in [41]. 

Uniqueness properties for local Dirichlet-to-Neumann maps associated with 
the wave equation are rather well understood (e.g., Belishev [4], Katchlov, Kurylev 
and Lassas [29], Kurylev and Lassas [30]) but stability for such operators is far 
from being apprehended. For instance, one may refer to Isakov and Sun E71 
where a local Dirichet-to-Neumann map yields a stability result in determining 
a coefficient in a subdomain. As for results involving a finite number of data in the 
Dirichlet-to-Neumann map, see Cheng andNakamura [14], Rakesh [34]. There are 
quite a few works on Dirichlet-to-Neumann maps, so our references are far from 
being complete: see also Cardoso and Mendoza |[T3l . Cheng and Yamamoto lfl5l . 
Eskin m-flS-lUa, Hech and Wang (23l, Rachele E3, Uhlmann fi_3] as related 
papers. 

The main goal of this paper is to study the stability of the inverse problem for 
the dynamical anisotropic wave equation. The approach that we develop is a dy- 
namical approach. Our inverse problem corresponds to a formulation with bound- 



ary measurements at infinitely many frequencies. On the other hand, the main 
methodology for formulations of inverse problems involving a measurement at a 
fixed frequency, is based on L 2 -weighted inequalites called Carleman estimates. 
For such applications of Carleman inequalities to inverse problems we refer for 
instance to Bellassoued [6], Isakov [26]. Most papers treat the determination of 
spatially varying functions by a single measurement. As for observability inequal- 
ities by means of Carleman estimates, see HI, (3, iflOl . 



Our proof is inspired by techniques used by Stefanov and Uhlmann 1411 . and 
Dos Santos Ferreira, Kenig, Salo and Uhlmann [17]. In the last reference, an 
uniqueness theorem for an inverse problem for an elliptic equation is proved fol- 
lowing ideas which in turn go back to the work of Calderon [ 121 . The heuristic 
underlying idea is that one can (at least formally) translate techniques used in solv- 
ing the elliptic equation 

3 t 2 + A g 

(which is the prototype of equations studied in lfT7l ) to the case of the wave equa- 
tion 

by changing t into it. Our problem turns out to be somehow easier because we 
don't need to construct complex geometrical solutions, but can rely on classical 
WKB solutions. 



1.1 Weak solutions of the wave equation 

Let (Ai , g) be a (smooth) compact Riemannian manifold with boundary of di- 
mension n > 2. We refer to [28] for the differential calculus of tensor fields 
on Riemannian manifolds. If we fix local coordinates x = [x\, . . . , x n ] and let 

■^r, . . . , gj- denote the corresponding tangent vector fields, the inner product 
and the norm on the tangent space T X A4 are given by 

n 

g(X,Y) = (X,Y) g = Y,m<xjPk, 



j,k=l 

n a n a 

\x\ g = (x,x)l/\ x = J2 ai —, r = X>^- 

i=l ' i=l 

If / is a C 1 function on M., we define the gradient of / as the vector field V g / 
such that 

X(f) = (V g f,X) R 



for all vector fields X on M. In local coordinates, we have 



df d 



^— ' dxi dx,; 



(1.4) 



«,i=i 



1 /9 

The metric tensor g induces the Riemannian volume dv™ = (det g) ' dxi A • • • A 
dx n . We denote by L 2 (M) the completion of C°°(M) with respect to the usual 
inner product 



</l,/2 



fi(x)f 2 (x)dv" f 1 J 2 GC 00 (M). 



M 



The Sobolev space H l (M) is the completion of C°°(M) with respect to the norm 

II " \\m(M)> 



m(M) ~ WJ \\L 2 (M) + II V /IIl 2 (>!) 



The normal derivative is given by 



d v u := V g u • v = 2_^ S 3 

j,k=i 



k d u 
3 dx k 



(1.5) 



where u is the unit outward vector field to dM. Moreover, using covariant deriva- 
tives (see |22]), it is possible to define coordinate invariant norms in H k (M), 
k>0. 

Let us consider the following initial boundary value problem for the wave equa- 
tion 

' (df- A s + q(x))v(t,x) = F{t,x) in (0,T) x M, 



v(0,x) = 0, d t v(0,x) = 

v(t,x) = 
The following result is well known (see |[24l ). 



in M, 

on (0,T) x dM. 



(1.6) 



Lemma 1.1 Let T > and q G L°°(M), suppose that F £ Jf , with Jtf = 
L 1 (0, T; L 2 (M)). The unique solution v of dl-6b satisfies 

veC 1 (0,T;L 2 (M))nC(0,T;H^(M)) 

and the mapping F \— > d v v is linear and continuous from ffl to L 2 ((0, T) x dM). 
Furthermore, there is a constant C > such that 

PAt, -)\\ L 2 {M) + ||Vu(t, -)\\ L 2 {M) < C \\F\\ L i {0yT . L 2 {M)) , (1.7) 

\\dv v \\L 2 ((0,T)xdM) <C\\F\\j£,. (1.8) 



A proof of the following lemma may be found for instance in ll3TTl . 

Lemma 1.2 Let f G H 1 ((0, T) x &M) &<? a function such that /(0, x) = 0/or a/Z 
x G <9A4. There exists an unique solution 

u£C 1 {0,T;L 2 (M))r\C{0,T;H 1 {M)) (1.9) 

fo ?/je problem rti.2D . Furthermore, the map f i— ► c^u jj linear and continuous 
from fl'HCO.T) x &M) mfo L 2 ({0,T) x 3A4). 

Therefore the Dirichlet-to-Neumann map A g)9 defined by (11.31 ) is continuous. We 
denote by ||A gj(? || its norm in £ (^((f^T) x 9A^); L 2 ((0,T) x &M)). Our last 
remark concerns the fact that when g is real valued, the Dirichlet-to-Neumann map 
is self-adjoint; more precisely, we have 

A* — A - 

This simple fact will be proven in section 12 We denote 

A g = A gi0 

the Dirichlet-to-Neumann map when there is no potential in the wave equation. 

1.2 Statement of the main results 

In this section we state the main stability results. Let us begin by introducing an 
admissible class of manifolds for which we can prove uniqueness and stability re- 
sults in our inverse problem. For this we need the notion of simple manifolds fiTTl . 

Definition 1 We say that the Riemannian manifold (M , g) (or more shortly that 
the metric g) is simple, if dM is strictly convex with respect to g, and for any 
x G M., the exponential map exp^, : exp~ 1 (A / l) — > Ai is a diffeomorphism. 

Note that if (M. , g) is simple, one can extend (M. , g) into another simple manifold 
Mi such that M d M\. 

Let us now introduce the admissible set of potentials q and the admissible set 
of conformal factors c. Let Mq > 0, k > 1 and e > be given. Set 

£(M Q ) = [q G H\M), \\q\\ HHM) < M } , (1.10) 

and 

V(M ,k,e) = 

{c G C°°{M), c > inM, ||1 - c\\ cl{M) < e, \\c\\ ck{M) < M } . (1.11) 
The main results of this paper can be stated as follows. 



Theorem 1 Let (M , g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2, let T > Di&m g (M), there exist constants C > and K\ £ 
(0, 1) such that for any real valued potentials q\, q 2 €E =S(Mo) such that q± = q2 
on the boundary dM , we have 

llffi -<&\\i?{M) ^ C|l A g,9i — A g,g 2 ir i (1-12) 

where C depends on M,T, Mq, n, and s. 

As a corollary of Theorem [Q we obtain the following uniqueness result. 

Corollary 1 Lef (M ,g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2, let T > Diam g (.M), let q\, q<i E J2(Mo) be real valued po- 
tentials such that q\ = q<± on dM. Then A gjgi = A gi g 2 implies q\ = q2 everywhere 
in M. 

Theorem 2 Let (M. , g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2, letT > Diam g (A / l), there exist k > 1, e > 0, < «2 < 1 
and C > such that for any c G ^(Mo, k, e) with c = 1 near the boundary dM., 
the following estimate holds true 

||l-c|| L2(A1) <C||A g -A cg |r 2 (1.13) 

where C depends on (M, g), Mq, n, e, k and s. 

As a corollary of Theorem |2l we obtain the following uniqueness result. 

Corollary 2 Let (M , g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2, let T > Diam g (A / l), there exist k > 1, e > 0, such that for 
any c £ ^{M^^k^e) with c = 1 near the boundary dM, we have A cg = A g 
implies c = 1 everywhere in M. 

1.3 Spectral inverse problem 

For q G J2(Mq) and g > 0, we denote by A q the unbounded operator A q = 
-A g + q with domain @(A q ) = H%(M) n # 2 (X). 

The spectrum of ^ consists of a sequence of eigenvalues, counted according 
to their multiplicities: 

< Ai,, < X 2 , q <..< Xk, q <■■■ 

with linifc^oo Xk t q = oo. The corresponding eigenfunctions are denoted by (4>k,q)- 
We may assume that this sequence forms an orfhonormal basis of L 2 (M). 

8 



In the sequel C denotes a generic positive constant depending only on M. and 
Mq (Mo is given by dl.lOft "). Since 4> k>q is the solution of the following boundary 
value problem 

(-Ag + q)<f> = X k>q (f> in M 
= 0, on 8M, 

classical H 2 {M. ) a priori estimates imply 

Hk, q \\ H , {M) < C\l[ 2 q U Kq \\ LHM) = C\l[ 2 q , * = 0, 1,2. (1.14) 

Therefore 

\\dv4>k,q\\HVZ(dM) — C^k,q- 

On the other hand, by Weyl's asymptotics, there exists a positive constant C > 1 
such that 

C- X k 2tn < \ Kq < Ck 2/n . (1.15) 

Here C can be chosen uniformly with respect to q provided < q{x) < M for 
x G M.. Therefore we have 

\\d u (t)k,q\\ H l/2( dM -} < Ck 

We fix r such that n/2 + l<r<n + l and it follows that 

(k~ 2r,n \\d v <t)k, q \\ H i/2 {dM) ) ei 1 . 

We recall that i 1 is the Banach space of real-valued sequences such that the cor- 
responding series is absolutely convergent. This space is equipped with its natural 
norm. 

Let uj = (oj k ) be the sequence given by ui k = k~ 2r ' n for each k > 1. We 
introduce the following Banach spaces 

{/i = (h k ) k ; h k G H 1 / 2 {8M) , k > 1, and (u k ||/»fc||ffi/3(aA<)) G i 1 } . 

and 

t (C) = {y = [y k )k] y k eC,k>l, and (u k \y k \) k G I 1 } . 

The natural norms on those spaces are 

V4ll{HV*(dM)) = Z2 uk W h k\\H 1 / 2 (8M) 
fe>l 



and 

ll y 'k(c) =£wfcll/fcl- 
fe>l 

We will apply Theorem [T]to prove the following result. 

Theorem 3 Let (M. , g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2. There exist C > and k% 6 (0, 1) such that the following 
estimate holds 

hi-q2\\ L 2 {M) <Ce^ (1.16) 

for any non-negative q\,q2 G £?(Mq) which are equal on the boundary dM, where 

6 = \\\ - ^52^1(0 + 11^^91 ~~ 9 u (j) q2 \\ e i ( H W*(dM)) 
is assumed to be small and d u (j) qj = ipv<j>k,qj) k > 3 = 1) 2. 

Theorem [3] is an extension of a result in |[T6l which is itself a variant of a theo- 
rem in El-||7). To the best of our knowledge, [2] is the first result in the literature 
concerned with stability estimates for multidimensional inverse spectral problems. 



The outline of the paper is as follows. In section 2 and 3 we collect some of the 
formulas needed in the paper. In section 4 we construct special geometrical optics 
solutions to the wave equation. In section 5 and 6, we establish stability estimates 
for related integrals over geodesies crossing Ai and prove our main results. In 
section 7 we prove Theorem 3. 

2 Preliminaries 

In this section we collect formulas needed in the rest of this paper. We denote 
by divX the divergence of a vector field X G H l (TAi) on M, i.e. in local 
coordinates, 

1 ^./ r^~ , .. ^ a 



v t=l i=\ 

If X G iJ 1 (TA^) we have the divergence formula 

/ divXdv"= I (X,v) da™" 1 (2.2) 

JM JdM 



10 



and for / G ^(M) Green's formula reads 

[ divXfdv^ = - [ (X,V g f) dv n g + [ (X,v)fda%- 1 . (2.3) 

JM JM JdM 

Then if / G H l (M) and w G H 2 (M), the following identity holds 

/ A gW /dv™ = - / (VgU/.Vg/) dv™+ / a^/da^ 1 . (2.4) 

Let /i, /2 G i? 1 ((0, T) x <9.M), we denote by u\, respectively by u 2 , the solutions 
to (11.21 ) with potential q and Dirichlet datum f\, respectively q and Dirichlet datum 
J2- By Green's formula, we have 

/ Wi£<K~ 1= / A g n 1 ^dv^+ / <V g «i,V^>_dvS 
•/&M ./jvp — v — ' jm 

=u\ qU2 

= / iti A g u 2 dv™ + / (V g Mi,V g it 2 )_ dv™ 



This shows that 



A4 JM 

hT~f 2 da n g - 1 . 
dM 



A* — A - 

g>9 ~~ §><?' 



In particular, this implies that A gi9 is selfadjoint when q is real-valued (and there- 
fore A g ). From now on, we will suppose the potential to be real-valued. 

For x G M. and G T X M we denote by ry x g the unique geodesic starting at 
the point x in the direction 9. We denote 



SM = {(x,6)eTM; \0\ g = l} , 
S*M = {(x,p)eT*M;\p\ g = l] 



the sphere bundle and co-sphere bundle of M.. The exponential map exp^, : 
T X M — > M is given by 

ex PxO) =lx,v{\ v \ g v ) =lx,v(rv), r = \v\ g . (2.5) 

A compact Riemannian manifold (j\A , g) with boundary is a convex non-trapping 
manifold, if it satisfies two conditions: 



11 



(a) the boundary dA4 is strictly convex, i.e. the second fundamental form of the 
boundary is positive definite at every boundary point, 

(b) for every point x G A4 and every vector 9 G T X M., 9 / 0, the maximal 
geodesic 7^0 (i) satisfying the initial conditions 

7x,0(O) = x and 7^0(0) = 9 

is defined on a finite segment [r_(x, 9), T+(x, 9)]. We recall that a geodesic 
7 : [a, b] — > M is maximal if it cannot be extended to a segment [a — e± , 6+ 
£2], where £j > and e\ + £2 > 0. 

The second condition is equivalent to all geodesies having finite length in M. . An 
important subclass of convex non-trapping manifold are simple manifolds. Recall 
that a compact Riemannian manifold (.M , g) which is simple satisfies the following 
properties 

(a) the boundary is strictly convex, 

(b) there are no conjugate points on any geodesic. 

A simple n- dimensional Riemannian manifold is diffeomorphic to a closed ball in 
W l , and any pair of points on the manifold can be joined by an unique minimizing 
geodesic. 

In the rest of this article, C will be a generic constant which might change 
from one line to another, but which only depends on the quantities allowed in the 
statement of the theorems (namely the quantities involved in the sets =2, ^, the 
manifold (Ai, g), the dimension n, the final time T and the Holder exponents Kj). 

3 The geodesical ray transform 

We introduce the submanifolds of inner and outer vectors of SM 

8±SM = {{x,9) e SM, x g dM, ±(9,v(x)) < 0} (3.1) 

where u is the unit outer normal to the boundary. Note that d+SA4 and d-SM. 
are compact manifolds with the same boundary S{dM), and 8SA4 = d+SM. U 
dSM. For (x,9) € d + SM, we denote by 7^ : [0,t + (x,9)\ — > M the 
maximal geodesic satisfying the initial conditions 7^0(0) = x and 7^61(0) = 9. 
Let C°°(d + SM) be the space of smooth functions on the manifold d + SA4. The 
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ray transform (also called geodesic X-ray transform) on a convex non-trapping 
manifold M is the linear operator 

X:C QO {M)^C 00 {d + SM) (3.2) 

defined by 

lf(x,0)= [ T f(rfx,e(t))db. (3.3) 

Jo 

The right-hand side of (13.31 ) is a smooth function on d + SM because the integration 

bound t + (x, 9) is a smooth function on d+SM, see Lemma 4.1.1 of IT381 . The ray 

transform on a convex non-trapping manifold M can be extended to a bounded 

operator 

Z:H k (M)^ H k {d + SM) (3.4) 

for every integer k > 1, see Theorem 4.2.1 of (38J. 

The Riemannian scalar product on T X M induces a volume form on S X M de- 
noted by duj x (9) and given by 

n 

doj x (9) = ^2(-l) k 9 k d6 1 A • • • A dO k A • • • A dd n . 
fe=i 

We introduce the volume form dvg™ -1 on the manifold SM. 

dvjf- 1 ^) = \&j x {P)/\ dv^| 

where dv" is the Riemannnian volume form on M.. By Liouville's theorem, the 
form dvg n_1 is preserved by the geodesic flow. The corresponding volume form 
on the boundary dSM = {(x, 9) £ SM, x £ dM] is given by 



r 2 g n - 2 = \du; x (9) A da^l 



where der™ 1 is the volume form of dM. . 

Let L?(d + SM) be the space of real valued square integrable functions with 
respect to the measure ii(x,9)da 2n ~ 2 with density /J,(x,9) = \(6,i/(x))\. This 
Hilbert space is endowed with the scalar product given by 



( u ' v )Ll(d + SM) = / u{x,9)v(x,9)n{x,9)da 2 

M Jd+SM 



r 2 g n - 2 . (3.5) 

Id+SM 
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The ray transform X is a bounded operator from L (M) into L (d + SM.) and its 
adjoint X* : L\ (d + SM ) — > I? (M ) is given by 



Ttj){x) = / if>*(x,9)du x (0) (3.6) 

JS X M 

where ifi* is the extension of the function ^ from d + SAi to 5.M constant on every 
orbit of the geodesic flow, i.e. 

r(x,e) = ^ x ,e(r+(x,e))). 

Let (Ai, g) be a simple metric, we assume, as we may, that {M., g) extends smoothly 
into a simple manifold such that M.\^> M.. Then there exist C\ > 0, C2 > such 
that 

Ci ||/|| i2(M) < \\X*X(f)\\ Hl{Mi) < C 2 \\f\\ L , {M) (3.7) 

for any / G L?{M.), see Theorem 3 in [40]. If V is an open set of the simple Rie- 
mannian manifold (.Mi, g), the normal operator X*X is an elliptic pseudodifferen- 
tial operator of order -lonF whose principal symbol is a multiple of |£| (see 
(32H401). Therefore there exists a constant C k > such that for all / e H k (V) 
compactly supported in V 

\\X*X(f)\\ Hk+1{Mi) <C k \\f\\ HHv) . (3.8) 

4 Geometrical optics solutions 

We will now construct geometrical optics solutions of the wave equation. We ex- 
tend the manifold (M,g) into a simple manifold M. 2 §> M and consider a simple 
manifold (7Wi,g) such that M. 2 3> M\- The potentials qi,q 2 may also be ex- 
tended to M. 2 and their ff 1 (7Wi) norms may be bounded by Mq. Since q\ and q 2 
coincide on the boundary, their extension outside Ai can be taken the same so that 
qi = q 2 inM 2 \M\. 

Let us assume for a moment that there exist a function ip G C 2 (Ai) which 
satisfies the eikonal equation 

\^K = t *££: = 1. Vx € A< 2 (4.1, 

and a function a G if 1 (IR, H 2 (M.)) which solves the transport equation 

R+t^^ + W)- - «*«."* (4.2) 

14 



with initial or final data 

a(t, x) = 0, Vx e M, and i < 0, or t > T. 
We also introduce the norm H-^ given by 

ll a ll* = \\ a \\H 1 (0,T-H 2 (M)) + \\ a \\H 3 (0,T;L 2 (M)) ■ 



(4.3) 



(4.4) 



Lemma 4.1 Le? g G L°°(M.),for any A > 0, the equation 

(<9 t 2 - A g + q{x))u = 0, in M T ■= (0,T) x M, 
u(k,x) = dtu(n,x)=0, K = 0,orT 

has a solution of the form 

u(t, x) = a{t, x)e iAW(a;) - t) + v x (t, x), (4.5) 

such that 

ueC 1 (0,T;L 2 (.M))nC(0,T;F 1 (.M)), (4.6) 

and where v\(t, x) satisfies 

v x (t, x) = 0, V(t, x) e (0, T)xdM, 

v\(k,x) = 0, dtV\(K,x)=0 x £ Ai, k = or T 
and 

a \\v x (t, -)WmM) + HMt, ^Wlhm) + ll v ^(t, OIL^m) < c Ml ■ ( 4 - 7 ) 

The constant C depends only on T and M (that is C does not depend on a and X). 
Proof . We set 

k(t,x) = -(d?-A g + q)(a(t,x)e iX ^- t) y (t,x) € (0,T) X M. (4.8) 
To prove our Lemma it would be enough to show that if v solves 
' (<9 2 - A g + q) v(t, x) = k(t, x) in (0,T) x M, 

v(k, x) = 0, d t v( K , x) = in M, r = 0, or T (4.9) 

k v{t,x) = on (0,T) x dM, 
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then the estimates (14.71) holds. We shall prove the estimate for k = 0, and the 
k = T case may be handled in a similar way. We have 

- k(t, x) = e^M"*) (d 2 - A g + q(x)) (a(t, x)) 

di\) da a 
T j - 

j,k=i 



I z — ' oxj oxk 2 



n 



+ \ 2 a(t,x)e^^ I 1 - £ g jfc |^ • (4-10) 

Taking into account (14.11 ) and (14.21) . the right-hand side of (14.101 ) becomes 

k(t,x) = -e^M-t) (a t 2 - A g + q) (a(t,x)) = -e^W^feo^x). (4.11) 
where fe G ^(O^T; L 2 (M)) and satisfies 

IIMl 2 ((0,T)xA4) + ||<9t/co|lL2(( 0jT ) xA 4) < C ||a||„ . 

Since the coefficient q does not depend on t, the function 

w\(t,x) = / v x (s,x)ds 
Jo 



solves the mixed hyperbolic problem (14.91 ) with right-hand side 

h(t,x) = J k(s,x)ds = ir I k {s,x)d s (e iX ^- s A ds. 

Integrating by parts with respect to s, we conclude that 

C 

H fc lllL 2 ((0,T)xA4) - y H a ll* • 

By Lemma [TTT1 we find 

v x GC 1 (0,T;L 2 (M))nC(0,T;H^(M)) (4.12) 

and 

C 
IM*, -)IIl2(A4) = \\ d tW\(t, -)llra(A4) ^ ^ H a ll* ■ ( 4 " 13 ) 

Since || k\\ L 2tr Q t)xM) — ^ ll a ll*' usrn g again the energy estimates for the problem 
(14.9b . we obtain 



\\dtv x (t,-)\\ L 2 {M) + \\Vv x (t,.)\\ L2{M) <C\\a\l. (4.14) 

The proof is complete. □ 
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Remark 1 In the construction of geometrical optics solutions, it is not necessary 
to assume that the potential is time independent. In the case where the potential q 
is also time dependent, one can proceed along the following lines. With the same 
notations, w\ satisfies the equation 

(Of - A g + q)w\ = k\ + / (q(t,x) - q(s,x))v\(s,x) ds. 

Jo 

If one uses Lemma 1.1 on the interval [0, r] one gets 

\\d t w\(T, -)\\vt(M) + l|V g t»A(r, -)IIl 2 (ai) 

C f T 

< -r \\ a L + CVf\\q\\ L °° / \\dtw(s, -)IIl2(a4) ds 

and Gronw all's inequality allows to conclude 

\\v\(t,-)\\l*(M) = II^A(r,-)||^(M) <jhL (l + Texp(CT 3 / 2 \\q\\ L ~)y 

We now proceed to construct a phase function tp solution to the eikonal equa- 
tion (I4.ll) and an amplitude function a solution to the transport equation (14.2I ). 



Let y £ dM.\. Denote points of J\A\ by (r, 9) where (r, 6) are polar normal 
coordinates in M.\ with center y. That is x = exp y (r6) where r > and 

e g s y Mi = {e 6 t v Mi, iei g = 1} • 

In these coordinates (which depend on the choice of y) the metric takes the form 

g(r,0) = dr 2 + g (r,6») 

where go(r, 0) is a smooth positive definite metric on S y M\. For any function u 
compactly supported in M., we set for r > and € S y M.\ 

u(r,9) = u(exp y (r9)) 

where we have extended u by outside M.. To solve the eikonal equation (I4.ll ) it 
is enough to take 

4>(x) = d g (x,y). (4.15) 

Then by the simplicity assumption, since y € .M2VM, we have ^ G C°°(A^) and 

^(r,0) = r = d g (x,y). (4.16) 
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We now proceed to the transport equation (14.21) . Recall that if f(r) is any function 
of the geodesic distance r, then 

Ag/(r) = / // (r) + ^-^/ / (r) (4.17) 

where a = a(r, 9) denotes the square of the volume element in geodesic polar 
coordinates. The transport equation (14.21 ) becomes 



Thus a satisfy 



da dip da 1_ _ 1 dadip 
dt dr dr 4 dr dr 



da da 1„ _ 1 da 
dt dr 4 dr 

Let G Q)°(R) and b G H 2 (d + SM), we choose a of the form 

a(t, r, 9) = a~ 1/4 0(t - r)b(y, 9). (4.20) 

A simple calculation shows that 

^(t,r,9) = a- 1 / 4 <P'(t-r)b(y,9). (4.21) 

and 

gf ( t , r , 0) = _ ^-5/4^^ (i _ r)6(y> 0) _ a -l/4^( t _ r ) 6 (y ; 0). (4.22) 

Finally, d4~22l and d472TT > yield 

9a, „. (9a. „. 1 ,_. „. 9a 

— (t, r, 9) + — (t, r, 9) = -- a - l a(t, r,9)—. (4.23) 

dt dr 4 ar 

If we assume that supp0 C (0, £o), with £o > small enough so that 

T > Diam g (7W) + 4e , (4.24) 

then for any x = exp (r#) G M., it is easy to see that a(t, r, 9) = if t < and 
t>T. 

Remark 2 If T > Diam g {M. ) + Ae$ and eg is e-close to g, then we also have 

T > Diam cg (.M) + 3e - 
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5 Stability estimate for the electric potential 

In this section, we complete the proof of Theorem [Q We are going to use the 
geometrical optics solutions constructed in the previous section; this will provide 
information on the geodesic ray transform of the difference of electric potentials. 

5.1 Preliminary estimates 

The main purpose of this section is to present a preliminary estimate, which relates 
the difference of the potentials to the Dirichlet-to-Neumann map. As before, we let 
<7i, <?2 £ .S(Mo) be real valued potentials. We set 

q = (qi-Q2)eH 1 (M). 

Recall that we have extended (71,(72 as H l (M. 2 ) in such a way that q = on 

M 2 \M. 

Lemma 5.1 There exists C > such that for any a\, a 2 G H 1 ^, H 2 (A4)) sat- 
isfying the transport equation (14.21) with initial data (14.31 ) the following estimate 
holds true: 

' / q(x)a\(t, x)a 2 (t,x) dv™d£ < 
Jm 

C(A- 1 + A 3 ||A g , (?1 -A g , g2 ||)||a 1 ||J|a2L (5-1) 

for any sufficiently large A > 0. 

Proof . First, if a 2 satisfies ( 14.21 ), (14.3b and A is sufficiently large, Lemma |4~T1 
guarantees the existence of the geometrical optics solutions u 2 

u 2 (t, x) = a 2 (t, z)e tt ^^ + V2,\(t, x), (5.2) 

to the equation with the electric potential q 2 

(<9 t 2 - A g + q 2 (x))u(t,x) = in(0,T) X M, 

u(0,-)=0, <9 t u(0,-) = inA^ 

where v 2 \ satisfies 

A IKaO, ■)\\ L 2 {m) + l|Vf 2 ,A(i, •)II L 2 ( _ M) < C ||a 2 |U 
V2,\(t,x) = 0, \/(t,x) G (0,T) x dM, 
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(5.3) 



and 

u 2 eC 1 (o,r ; L 2 (^))nC(o,T;ff 1 (M)). 

Let us denote by f\ the function 

f x (t,x) = a 2 {t,x)e iXWx) - t \ t G (0,T), x G AM. 
Let i> denote the solution of the following initial boundary value problem 
f (d?-A g + qi )v = 0, (t,x)E(0,T)xM, 

v(0, x) = 0, d t v(0, x) = 0, xeM, 

k u(t,z) = «2(t,ar) := / A (i,x), (i,x) G (0,T) x 8M. 
Taking w = v — u 2 , one gets 

' (d?-A g + qi (x))w(t,x) = q(x)u 2 (t,x) (t,x)e(0,T)xM, 



(5.4) 



iu(0,x)=0, d t w(Q,x)=0, 



x £ M, 



w(t,x) = 0, (t,x) G (0,T) x dAL 

Since q{x)u 2 G L 1 (0, T; L 2 {M)) by Lemma ITTTl we deduce that 

weC 1 (0,T;L 2 (M))nC(0,T;H 1 (M)). 

Therefore, we have constructed a particular solution u\ G C 1 (0, T; L 2 (A4)) n 
C(0, T; i? 1 (A4 )) to the backward wave equation 

{d 2 t - A g + q x (x)) ui (t, x) = 0, (t,x)€ (0, T) x M , 

u 1 (r,ar) = 0, <9 t ni(T,x) = x G X, 

having the form 

Ul (t,x) = ai {t,x)e iX ^^-^+vi, x (t,x), (5.5) 

corresponding to the electric potential q\, with 

A IKa(*, OIIl^) + HVvi.ACt, 0IL2 (M) < C ||oiL • (5-6) 

Integrating by parts and using Green's formula (I2.4I ). we find 

/ / (<9f-A g + gi(x))«nZidv?dt 
JO J M 

- I I d v wui da"" 1 dt. (5.7) 

Vo J9A4 



/0 JA1 



l(x)it2Wi dv" dt 
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Combining (15.71) with (15.4I ). we deduce 



nqu 2 ui dv™ dt 
.A 



r-T 
'0 JdM 



{K^ qi -K &q2 ){f x )g x da n g - l dt (5.8) 



where 



g x (t,x) = ai(t,x)e iA W a ')- t ), (t,x) G (0,T) x 9M 

It follows from (5M . d53T > and ( 1531 ) that 

" T , 

](x)a 2 (t, x)ai(t, x) dvg (it 

(A g)9l -A gi g 2 )(f x )(t,x)g x (t,x)da^- 1 dt 



JM 



'0 JdM 

(■T 



/0 7M 
l-T 



l{x)a 2 {t,x)v liX {t,x)e iX ^ dv"dt 
l{x)v 2X (t, x)ai(t, x)e~ iX ^-^ dv™ dt 



/0 J7K 



/0 J.M 

In view of (15.61 ) and (15.31 ). we have 

f f qa 2 vi,xe^-^ dv™ 
o Jm 



l(x)v 2 ,\(t, x)vi tX (t, x) dv™ dt. 



(5.9) 



<C \\a 2 (t,-)\\ L2{M) |K,a(V) 
< CA -1 H^IL H a ilL • 



lL 2 (A4) 



dt 



(5.10) 



Similarly, we deduce 

/ / q(x)a 1 (t,x)v 2 , x (t,x)e- iX ^- t ^dvUt 
'0 J M 

Moreover we have 

/ / q(x)v 2 , x (t,x)v 1>x (t,x)dvi l dt 

'o Jm 

On the other hand, by the trace theorem, we find 

cT 

_n-l 



< CA llailL II ct2 1 



< CA ||oi|| # ||a 2 |U • 



JdM 



( A g,gi- A g,9 2 )(/A)9Ad< Mt 



< ||A„. m -A 



g,9i A g,92ll H/A||_ffi((o,T)xaX) llfl , A||L 2 ((0,T)xaA4) 



< CA 3 ||oi||. ||a 2 |L IIA 

— II J-II^H ^11*11 



* ll J1 g,gi X1 g,?2 



K 



(5.11) 
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Inequality <|5J} follows easily from BM . (f57TOT) . <f57Tb . (ETTT) and (TSTTTb . This 
completes the proof of the Lemma. □ 



Lemma 5.2 There exist C > 0, j3 G (0, 1) smc/z that for any b G H 2 (d + SM\) 
the following estimate 



I f 

JSyMlJO 



■{y,o) 



u«/: -A gi92 || P HKy,-)|| H 2 (5 + A/ll) (5.12) 



(s,9)b(y,9)n(y,9)dsdoo y (9) 
<C||A 
holds for any y G dAi\. 
Here we used the notation 

S+Mi = {0 G V^i : M> < o}. 

We recall that fi denotes the density — (0, u(y)) g . 

Proof . Following (14.20b . we take two solutions to (14.21 ) and (14.31 ) of the form 

a 1 (t,r,9) = a- 1/ U(t-r)b(y,6), 
a 2 (t,r,9) = a- 1/ U(t-r)v(y,8). 

Now we change variables in (15.11 ). x = exp y (r9), r > and 6 G S y M.\, we have 



r 



/(x)ai(t, x)a2(i, x) dv™ dt 



o JM 



T+(y,o) 



'0 JS y XiJ0 
,T r rr+iyfi) 



<(r, B)ai (t, r, 9)a 2 (t, r, 0)a 1/2 dr dw y (0) dt 



(r, 0)</> 2 (i - r)6(j/, 0)/z(y, 0) dr d^(0) dt. 



'0 JSyMiJO 
By virtue of Lemma l5.ll we conclude that 



JSyMiJO 



-{y,0) 



q(r, 9)4? (i - r)6(y, 0V(y, 0) dr dw y (0) di 



<C(A _1 + A 3 ||A 



"•'/' A g! q 2 | 



h*(P)\\Kv,-)\\h*(s+Mi)- (5.13) 
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Since 0(t) = for t < or t > T, we get 

Jo JSyMiJo 



4> 2 (t)dt x 
oo / JSyMiJO 



(r, 0)0' (t - r)b{y, 9)fi(y, 9) dr duj y (9) dt 

r+(yfi) 



j(r,e)b(y,eMy,9)drdojy(e). (5.14) 



Combining (15.131 ) and ( 15.141 ), it follows that 
-(31,0) 



[ f T q(s,9)b(y,9)fx(y,9)dsduj y (9) 

JSyMiJO 

<c(2- + A 3 ||A g , gi -A g , g2 ||)||% 



''■>\\H*{S+Mi) ■ 



Finally, minimizing in A we obtain 



f r +V ' q(s,9)b(y,9My,9)dsdu; y (9) 

JSyMiJO 

<C\\A g>qi -A g>q2 f\\b(y,-)\\ H2{s + Mi) . 
This completes the proof of the lemma. □ 



5.2 End of the proof of the stability estimate 

Let us now complete the proof of the stability estimate in Theorem [TJ Using 
Lemma I5T21 for any y G dM\ and b G H 2 {d + SM) we have 



/ l(q)(y,9)b(y,9My,9)du; y (9) 

JSyMl 



< C ||A gigi - Ag.gall" ||6(y, •)llir»(s+A<0 
Integrating with respect to y G 9.M i we obtain 



/ l(q)(y,9)b(y,9)(9,u(y)) 

Jd+SMi 



daf- z (y,9) 



s 



< C||A. g)9l - A gj(?2 | 



l^ 2 (9+SXi) • ( 5 - 15 ) 
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Now we choose 

b(y,d)=l(l*l(q))(y,9). 

Taking into account (13.81 ) and (13.41) . we obtain 

ii^^qoiil^.mi) - cniAg>«i - A g,g 2 ii ikiifl-i(M) • 

By interpolation, it follows that 

< C \\Z*Z{q)\\ L 2 {Ml) \\q\\ mm 

<C\\Tl{q)\\ L2{Mi) 

<C||A g , gi -A g , 92 f/ 2 . (5.16) 

Using (13.71 ). we deduce that 

n n 2 s n n a a n/ 9 / 2 

II^IIl 2 (A4) - ° ll A g,qi ~~ A g,Q2ll 

This completes the proof of Theorem [Q 

Remark 3 In the proof of Theorem [Q we have used the time independence of the 
potential at two stages: 

1 . In the construction of the remainder term v\ (t, x) in the proof of Lemma l4TT1 
But as was noted in Remark this restriction may be bypassed. 

2. In equation (15.141) to get rid of the 4> 2 (t—r) term and obtain the ray transform. 
The adaptation to the time dependent case does not seem to be straightfor- 
ward. 

We leave the case of a time dependent potential as an open problem. 

6 Stability estimate for the conformal factor 

We shall use the following notations. Let c G ^{Mq, k, e), we denote 

Q (X) = 1 - C(X), Ql (x) = C n / 2 (x)-l, Q 2 ( X ) = C n / 2 - 1 (x)-l, 

g(x) = Ql (x) - q 2 (x) = c n l 2 ~\x) (c(x) - 1) . (6.1) 

Then the following holds 

\\qj\\ c i(m) - c WqoWchm) > i = 1 > 2 - ( 6 - 2 ) 

The first step in our analysis is the following result. 
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Lemma 6.1 Let c G C°°(A4) be such that c = 1 near the boundary dAi. Let u\, 
U2 solve the following problem in (0, T) x Ai with some T > 



f (<9 t 2 - A g )m = 0, in (0, T) x M, 



Ul (0,-) =3 t -u(0,-) =0 inM, 

I ui = /i, on (0,T) x &M, 

' (3f - A cg )u 2 = 0, m (0,T)x7W, 



u 2 (T,-) = d t u(T,-) =0 /« A4, 

I u 2 = / 2 , on (0,T) x9^ 

where /& £ i/ 1 ((0, T) x (9.M), k = 1,2. Then the following identity 



(6.4) 



(6.5) 



(A g — A C g) ji j 2 diTg at = / ^ iV ~ .,.-/.., ,,, .._/m, ,, 
/o ^9A4 jo iM 



Q 1 (x)d t uid t u 2 dv™ dt 



52(x) (V g «i(t, x), Vg« 2 (t, x))_ dv"dt (6.6) 



JjVl 



g s 



holds true for any fj £ H 1 ((0, T) x dM), j = 1, 2. 

Proof . We multiply both hand sides of the first equation (16.41 ) by «2 ; integrating 
by parts in time and using Green's formula (I2.3I >- (I2.4I ) in (Ai, g), we obtain 

= f j (dfut - A g ui) u 2 dv"dt 
Jo Jm 

nd t uid t U2 dv™ g dt+ / qi (x)d t u 1 d t U2 dv™ di 
m "Jo Jm 



+ 



JM 






v'fe 



#»1 *fc\ dy n dt 

dxj dxk ' cg 



1 '" ■ " du\ du2 



£ 



T i fe . 



r 



02 W ,/.«,« o 



dv"dt- / / S^ui/a da? -i dt. 

Jo JdM 
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Another integration by parts in time and an application of Green's formula in 
(M , eg) yield 

= f [ «i (d?u 2 - A cg u 2 ) dv" dt 
Jo Jm 

+ 11 d u u 2 fi da™" 1 dt- i [ d u uj 2 da^dt. (6.7) 

JO JdM JO JdM 

Taking into account the facts that A cg is self-adjoint, that c = 1 on dJ\A and 
(dfu 2 - A cg u 2 ) = in (0, T) x M, it follows that 

/ / (A g -A qj )/ 1 / 2 d^- 1 dt = 

JO JdM 

(6.8) 
This completes the proof of the Lemma. □ 



6.1 Modified geometrical optics solutions 

As in the case of potentials, we extend the manifold (JvL, g) into a simple mani- 
fold M 2 =) M so that M 2 D M± 2> M with (M\,g) simple. We extend the 
conformal factor c by 1 outside the manifold Jvt; its C k {M.\) norms may also be 
bounded by Mq. Let ipi, tp 2 be two phase functions solving the eikonal equation 
with respect respectively to the metrics g and eg. 



IV7 / |2 V^ jk u <Pl ul Pl 1 |V7 / |2 V^ jk u W2VW2 r 

i v g H = E ^ w s w k = !» i v ^i cg = E ^ &^ = L 

j,k=i J iifc=i 

(6-9) 



Let a 2 solve the transport equation inR x M with respect the metric g (as given 
in section 4) 

j,k=l J 
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Let as solve the following transport equation in R x M. with respect to the metric 

eg 

j,k=l ■> 

= a 2 (t,x)ip (x,X), (6.11) 

and be such that 

||a 3 |U < C'eA 2 ||a 2 |U • ( 6 - 12 ) 

Let us explain the construction of a solution 0,3 satisfying (16.111 ) and (16.121) . To 
solve the transport equation ( 16.111 ) and ( 16.121 ) it is enough to take, in geodesic polar 
coordinates (r, 9) (with respect to the metric eg) 

a 3 (i,r,0;A)=cC 1/4 (r,0) / a 1 J\s,6)a 2 (s-r + t,s,e)w(s,6,\)ds, (6.13) 

Jo 

where a cg (r, 9) denotes the square of the volume element in geodesic polar coordi- 
nates with respect to the metric eg. Using that \\ifo(-, -M llc 2 OW) — CeX 2 and ( 16.131 ) 
we obtain (16.121) . 

Lemma 6.2 Let c G ^(Mq, k, e) be such that c = 1 near the boundary dM. Then 
the equation 

(<9 t 2 - A cg ) u = 0, in (0,T)xM, «(0, x) = d t u(0, x) = (6. 14) 

has a solution of the form 

u 2 (t,x) = \a 2 (t,x)e a ^W-V +a 3 (t,x;\)e iX ^ x ^ +v 2 , x (t,x) (6.15) 
A 

when A is large enough, which satisfies 

A ||V2,a(* } -)llz,2(7U) + ||VgV2,A(<, •)II L 2 ( _ M) + ||^«2,a(*, OHi^M) 

<C(eA 2 + A- 1 )||o 2 |U. (6.16) 

TTie constant C depends only on T and Ai (that is C does not depend on a, A and 

e). 

Proof . We set 

k(t,x) = -{d?-A cg ) (ja 2 (t,x)e lX ^-^+a 3 (t,x,X)e lX ^-A. (6.17) 
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To prove our Lemma it would be enough to show that if v solves 

(df - A cg ) v = k(t, x) (6.18) 

with initial and boundary conditions 

v(0,x) = d t v(0,x) =0, inM, and v{t,x) = on (0,T) x dM (6.19) 
then the estimates (16.161 ) holds. We have 



1 



k(t,x) = -e^- t >{d?-A cg ) 



A 



0-2 



n 



j,k=l 



jk d2pida2_ ^ (H 
] Xj dxk 



dxj dxk 2 



\ ^i 9xjdx k 

+ e lX ^- t] id 2 , - A cg ) a 3 



j,k=l 



ifc dtj)2 da 3 03 
dxj dxk 2 



a 



+ AW^- i Ml-Ecg^||g . (6.20) 

Taking into account (16.91 ) and (16. 10b . the right-hand side of (16.201 ) becomes 

-k(t,x) = I e ^-*) {d 2 t - A cg ) a 2 

+ 2ie iX ^~^ ((c- 1 - 1) <V g Vi, V g a 2 ) g + X -a 2 (A cg ^i - A g ^ 



j,k=i 



dxj dxk 2 



+ V*-*)(l-c- 1 )) 
+ e^ 2 -*) (^ - A cg ) a 3 . (6.21) 
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By (16.111) we get 

-k(t,x) = je lX ^-^{d?-A cg )a 2 

+ 2ie lX ^-V (V 1 - 1) (V g Vi, V g a 2 ) g + ^a 2 (A cg ^i - A g ^) ) 



IgiA^i-t)^ + giA^i-At)^ + e iKi>2-t) k2 _ 



A 



Since fej G L^O, T; L 2 (.M)), by Lemma [TTT1 we deduce that 

^GC 1 (0,T;L 2 (X))nC(0,T;F 1 (X)) 
and 

IK (V) II 1,2 (M) 



(6.22) 
(6.23) 



C 

< — 
~ A 



+ 



j \\ko{s, -)\\ L 2 {M) + ||fci(a, -)IIl2 ( a4) + IMs, -)IIl2(m) 

T ||ftfco(s, -)IIl2(A4) + H^lO, -)lll,2(M) + ll$M*> Olll^M) 



d,9 



C ( \ 

- X 1 A "° 2 "* + e ll a2 "* + eA2 H a2 "* 



1 



< C I eA + ~2 ] Ik/;? 
Moreover, we have 



I^IIl 2 ((o,t)xm) - c (t "° 2 "* + e " a2 ll* + eA2 " a2 l 



(6.24) 



(6.25) 



and by using again the energy estimates for the problem (16.18I )- (I6.19I ), we obtain 

I|Vva(<, OIL^m) + W d Mt, ')\\v(M) ^ c ( £A2 + \ 

This ends the proof of Lemma l6\2l □ 



(6.26) 



Lemma 6.3 There exists C > such that for any a\, a 2 G i? 1 (M, H 2 (Ai)) satis- 
fying the transport equation ( I6.J0I ) with rt4.JP the following estimate holds true 



J M 



Q{x)(a\a<i)(t,x) dvgdt 



^ ll^ollc(M) ( A 1 + eX3 ) IKL IMI* 

+ A 3 ||ai||J|a 2 ||J|A g -A cg || (6.27) 



for any sufficiently large A. 
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Proof . Following Lemma l6T2~l let u 2 be a solution to the problem (df — A cg )u = 
of the form 

u 2 (t,x) = \a 2 (t,x)e- iX ^- t) +a 3 (t,x;X)e- iX ^- t) + v 2 ,x(t,x) 
A 

where v 2jX satisfies (16.161) and a 3 satisfies (16.121) . Thanks to Lemma l4~T1 let u\ be a 
solution to the {df — A g )u = of the form 

ux(t, x) = ax(t, x)e iA ^ 1 -* ) + v 1>x (t, x), 
where v\ X satisfies ( I4.4I ). Then 



d t u 2 (t,x) = \d t a 2 (t,x)e- iX ^'^ + ia 2 (t,x)e~ iX{ ^- t) 
X 

+ d t a 3 (t, x; \) e - iX ^-V + iXa 3 {t, x, X)e~ iX ^-^ + d t v 2 , x (t, x) 
d tUl {t,x) = 9 t oi(t,a;)e u ^ 1 -* ) -iXa 1 (t,x)e iX ^ 1 - t ') + d t v ljX . (6.28) 



Let us compute the first term in the right hand side of (16.6I ). We have 

nQid t uid t u 2 dVg dt = A / / £iaia 2 dVgdi 
.A JO JM 

+ T Qi (d t aid t a 2 ) dv^dt-i / / Qiaid t a 2 dv" dt 

A Jo Jm ' Jo Jm 

+ \ I I Qid t a 2 d t v hX e- lX ^-^dwldt + i [ [ 6l (8^02) dv" dt 
A Jo Jm Jo Jm 

+ i f I Qia 2 d t v 1 , x e- iX ^-^dv g l dt+ [ [ ei dta 3 d t a 1 e iX ^ 1 -^ dv n g dt 
Jo Jm Jo Jm 

-iX f [ Q 1 dta 3 a 1 e lX ^ 1 -^dv2dt+ [ [ gi d tVl A a t a 3 e" iAW2 " f) dv"d£ 



Jo jm ' Jo Jm 

+ iX f f e^sdt^e^ 1 -^ dv"di + A 2 / f g^iasJ^-^dv^dt 
Jo Jm Jo Jm 

+ iX f f 01 dtv hX a 3 e- lX ^^dv^dt+ [ [ gxd t a 1 d t v 2t xe iX ^- t) dv^dt 

Jo Jm Jo Jm 

-iX i i Ql a 1 d t v 2 . x e iX ^-^dv n g dt+ [ I Qid t v liX dtV 2 , X dv2dt. (6.29) 
Jo Jm Jo Jm 

Thus, we have from (16.12I ). (16.161 ) and (14.41) the following identity 

nQi(x)d t uid t u 2 dv g dt = X / / Qi{x)(aia 2 )(t, x) dv g dt + Ji(A,e), 
A Jo Jm 

(6.30) 
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where 

\Jl(\e)\ < \\go\\c(M) I 1 + eA4 ) IM* ll°ill* • ( 6 - 31 ) 

On the other hand, we have 

V g ui = VgOie**^ 1 -^ + lAVg^aie^ - * 5 + V g i» ljA 

V g n 2 = \y g a 2 e- iX ^^ - i^VgV'ie"**^ 1 "* 5 
A 

- tAasVg^e-^^-*) + Vgase"^ 2 ^ + V g T7 2 , A . (6.32) 
and the second term in the right side of (16.61 ) becomes 

ng 2 (x) (V s ui(t,x),V g u 2 (t,x)} dv^dt 
M 

A 



with 



/ / Q 2 (x)(a 1 a 2 )(t,x)dv g i dt + J 2 (\,e) + J 3 (\,e) (6.33) 

JO JM 



M\,e)= + ~J J e2(V g ai,V g 02) g dvJdt 

-ill 62a 2 (V g a 1 ,V g ipi(x)) „ dv g di 

Jo J,M 

!/ / £> 2 ai(Vg0 2 ,VgV'i) g dv^dt 



+ i 



and 



+ T /7 ^e- iA( ^- t) (V g a 2 ,V g? ; 1)A ) g dv^dt 

A Jo JM 

-iff ^a 2 e- iA ^- 4 )(V gVl , A ,V g ^i) K dv^dt 
JO JA4 



M\,e)=-i\[ f Q2ase iX ^-^(V g ai,VM R dv2dt 
Jo Jm 

+ J J 2 e a ^-« (V g0l , V g a 3 ) g dv™ d* 
+ / / ^e aWl - t) (Vgai,VgtJ 2!A ) K dv^dt 

JO JA4 ■ ■ b 

+ A 2 f f ^a 1 a 3 e^ 1 -^ ) (V g Vi,V g ^ 2 ) g dv^dt 
Jo JM 
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Jo Jm 

+ iX f f ^aie aWl-t) (V g Vi,V g V2} B dv"dt 

Jo Jm 6 6 

rp 

-ixf f Q2a 3 e- iX ^- t Uv g v 1A ,V 9 ij2) R dv n g dt 
Jo Jm • • s 

rp 

+ 11 Q2e- lX ^~ t) (V g vi,x,V g a 3 ) g dv n g dt 
Jo Jm 

r-T 



+ 



/ / Q2{V g v 1:X ,V g v 2 ,x) s dv"di. 

JO JA4 



'0 JA4 
From (16. 12b . (16.161 ) and (l4~4l ). we have 

\M\e)\ + \M\,e)\ < ||eo|| c(M) (l + eA 4 )||a 2 ||J|ai| 

Taking into account (16.6I ). (16.301) and (16.331) . we deduce that 

/ / (A g -A cg )f 1 J 2 da-- 1 dt 
Jo JdM 



(6.34) 



= A/ / Q{x){a 1 a 2 )(t,x)dvyt + J 1 (\,e) + J 2 (\,E) + J 3 (\,e). (6.35) 
Jo Jm 

In view of (16.341 ) and (16.311) . we obtain 



T , 



g(x)(aia 2 )(t,x) dv™dt 



o JM 



< \\Qo\ 



C(M) 



A 1 +eA ,i ) I oil \\a 2 \ 



This completes the proof. 



+ A 3 ||ai||J|a 2 ||J|A g -A C g||. (6.36) 



D 



6.2 Stability estimate of the geodesic ray transform 

Lemma 6.4 Let Mq > 0. There exist C > and /3j > 0, j = 1, 2, 3, such that for 
any b G H 2 (d + Sj\A.i) the following estimate 



I 



d+SM 



T{ Q ){y,9)b{y,e){e,u{y)) dal n -\y,6) 



< ((A-^ + ^ 2 )||eo|| c i ( M)+ A/33 l|A g -A cg ||)||6|| H2(a+5Mi) . (6.37) 
holds for any A Zarge. 



32 



Proof . Following (14.201) . we take two solutions of the form 

a 1 (t,r,e) = a- 1 /U(t-r)b(y,6), 
a 2 (t,r,9) = a- 1 /U(t-r)fi(y,6). 
Now we change variable in ( 16.271 ). x = exp y (r0), r > and 9 £ S y M\. Then 



T 



g(x)ai(t, x)a,2(t, x) dv™ dt 



o JM 




r+{yfi) 



'0 JSyMiJO 
rT r rr+(y,6 



g(r, 9)ai (t, r, 0)a 2 (t, r, 9)a 1/2 dr du y (6) dt 



g(r, 9)cf> 2 (t - r)b(y, 6)n(y, 9) dr du y {9) dt. 



'0 JSyMiJO 
We conclude that 



JSyMiJO 



r+(y,0) 



g(r, e^it - r)b(y, 9)fi(y, 9) dr du y (9) dt 

i \ 3 || A 

C(M) 



< C[ (A l + eA 3 ) ||eolU A ^ + A 3 ||A,.«i - A. 



,91 J1 g, 92 I 



H 3 



(R)\\Kv>')\\h*(S+Mi)- ( - 63 ^ 



Since 4>(t) = for t < or t > T, we get 

'■r+{yfi) 



oc 



J5 s MiJ0 

OO 

— oo 



e(r, 0)<^(t - r)b(y, 6)fi(y, 9) dr doo y (9) dt 



r+(y,9) 



cf> 2 (t)dt)x / Q(r,0)b(y,0)ii(y,9)dr6u v (p). (6.39) 

/ JSyMiJO 

Combining (16.391 ) and (16.381) . it follows that 
-(31,0) 

g(s,9)b(y,9My,9)dsduj y (9) 

< C^X- 1 + eA 3 ) ||0o|| c(M ) + A 3 ||A g)91 - A g>92 || J ||%, 011^(5+^0 
Integrating with respect to y £ dM. i we obtain 



SyMlJO 



&< n ~\y,e) 



f l(g)(y,e)b(y,9)(e,v(y)) 

Jd+SMi 

< ( (A -A + eA &) \\ eo \\ cl(M) + A^ || Ag - A cg || ) \\b\\ H2 



(d+SMi) 



(6.40) 
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This completes the proof of the lemma. □ 

6.3 End of the proof of Theorem |2] 

This subsection is devoted to the end of the proof of Theorem |2j We need the 
following known result (see [41] proposition 4.1). 

Lemma 6.5 Let c £ C°°(M) be such that ||1 — c|| C (^\ < e. Then there exists 
C > such that 

\\d g — dcg\\c(dMxdM) — C l|A g — A cg || , (6.41) 

with some < p, < 1 depending only on the dimension n. 

From this lemma, we can derive the following estimate 

Corollary 3 There exists a constant C > such that the following estimate holds 
true 

\\T*l{ Q )\\%t {Mi) < C (\\e\\ 2 eHM) + ||A g - A cg f) (6.42) 

with some < p, < 1 depending on n only. 
Proof . Linearizing near g, we get, as in iBTI 

d g (x,y) -d cg (x,y) = -l(g)(x,y) +K(g)(x,y), \/x,y G dM, 

where, with some abuse of notation, l(g)(x,y) stands for I(g)(x,6) with 6 = 
exp~ 1 (y)/ |exp~ 1 (y)|. The remainder term lZ(g)(x,y) is nonlinear and satisfies 
the estimate (see [21]) 

\K{o)(x,y)\ < Cd g (x,y) \\q\\c^m) > Vx '^ G dM > 
with C > uniform incifO<e^Cl.By Lemma [631 we have 

\l{Q){x,y)\ < C (d g (x,y) \\g\\ 2 c i (M) + ||A g - A cg f) , Vx,y e dM. 

Apply 1* to both sides, and use the estimate ||2^*(/)||^cx)(_ A/1 ) < C WHIl^IM ) io 
get 

11^(^)11^(^0 < O (llellci(M) + HAg - A cg r) . (6.43) 

Since g vanishes outside M with all derivatives and T*X is a pseudodifferential 
operator of order —1, we have 

\\Z 2(0)1111™+! (Mi) — Cm \\q\\h™(M) 
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for all integers m. Using interpolation, we get 

1 2/3 n—*—/ xnl/3 



2/3 
L°°(Mi) 



<C||X*X(^)||^_ V (6.44) 



Therefore, (16.441 ) and ( 16.43b imply 

||X*X(^)||5/ 2 2 (Mi) < C (|M& (M) + ||A g - A cg r N 
This completes the proof. D 

Let us now prove Theorem [2] We choose 

b(y,8)=l(ri(q))(y,e) 
and obtain 

\\Z*Z(Q)\\l HMl ) 

< C ((A-ft + eA*) [|eo[| C i (A<) + A* ||A g - A cg || ) \\T*1(q)\\ hHMi) . 

By interpolation we have 

\\l*l(6)\\ 2 H HMl ) 

<c\\ri(6)\\ L2{Mi) \\ri(6)\\ H2{Mi) 

< c((A-/fc +£A ^) [|eo|| cl(A4) + a^ || Ag - A cg || ) 1/2 \\i*i( e )\\f 2{Mi) . 

Using (16.421 ) we obtain 

\\Q\\h(M) < C{(^ [ + *>&) lleoll^M) + ^ l|A g - A cg f ) . 
Since 

II II ^ /-Ml II s- n\\ ll 4 /5 II II 1/5 //on 1 1 4/5 

\\eo\\ci(M) s o neoiii/~/2+i+€ (-M ) s o lieolli3(^() I l^o IIh«(m) - u H^oIIl 2 (a4) 

we obtain 

\\Q\\ 2 L 2 iM) < c({\-& + e\&) \\eo\\h {M) + A^ ||A g - A cg | 
Minimising (A _/?1 + e\^ 2 ) with respect to A > 0, we get 

C/ii ii2 ^ \\ \\1 ^>- s~i ( y 1 1 1 1 2 | •'-f 1 1 a a 1 1 U 

\\Q0\\L 2 (M) - \\&\\l 2 (M) - ° \ £ \\80\\l2(M) + °e H A g ~ A cgll 

for e > small enough we conclude and obtain dl - 13b - 
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7 Proof of Theorem g] 

Let q G L°°(A4), we first define the elliptic Dirichlet-to-Neumann map; let <r(A q ) = 
{^k,q} be the spectrum of A q and p(A q ) = C \ cr(A q ) be the resolvent set of A q . 
From well known results (e.g., ||3"T1 ), for any z G /o(Aj) and /i G H 3 / 2 {dM.), the 
nonhomogeneous boundary value problem 

(— A g + q)u = zu, inM 
u = h, on dM. 

has an unique solution in u q> h G H 2 {M.) and the Dirichlet-to-Neumann map 

n g , g (z) : / -)■ 0z/%/i|a.M 

defines a bounded operator from H 3 / 2 (dM) to i? 1 / 2 (3A / f ). We fix T > Diam g .A/f 
and consider the following function space 

«l = {/ £ H 2n+A (0,T;Hl(dM)); d(f(0,-) = 0, < j < 2n + 3 } , 
and the operator 

*w =et^2^« r siny ^" s) (-^ (n+2) /(^),^, fc )d,, 

fc>l A ?,fe - 70 V A 9. fc 

where (-,-) denotes the L 2 (<9.M)-scalar product. Then M gA defines a bounded 
operator from Hi to H 2 = L 2 {0, T; H s {dM)). 

We will need in the sequel the following three lemmas. Their proof can be 
found in [2] or can be deduced easily from the results in this reference (see also 
(HI). We fix < s < i. 

Lemma 7.1 Let q G L°°{M). Then for any m > §, h G H 3 / 2 (dM) and z G 
p(A g ), we have 

— n M (z)/i = -m! ^ - z )m+l ( h > d ^k, q )du<t)k, q . 

Lemma 7.2 Let N be a non negative integer and let q\, q<i G L°°(M) satisfy 
< <7i , 92 < ^o / or some positive constant M. Then there exists a positive 
constant C, depending only on Ai and Mq, such that 



d j 

n g,<?i(^)- n g,g 2 (^) 



- I , i-2, » -2 < ° aw ^ < j < iV, 



ur * 



03* 

w/iere || • [| denotes the norm in C(H 3 / 2 (dM.); H s (dM)). 
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Lemma 7.3 For each f £ Hi, we have 

n+l 

A Lf = E 



3=0 






U g,q( Z ) 



|z=0 



-d 2 t f)+M z J, 



(7.1) 



where Ag i? z's the restriction ofA g ^ q to ~Hi. 



First, we remark that for q G L°°{M) and / G %i the problem (11.21 ) has an unique 
solution 

€L 2 (0,T;# 2 (Al))n# 2 (0,T;L 2 (.M)) . (7.2) 



U 



Moreover A gj(? is a linear and continuous map from %\ into %• Indeed, for / G 
%i, let v solve the problem 



Then 
Furthermore 



f (a?-A g )t; = 0, in (0,T)xM, 

v(0,-) = 0, a t v(0, •) = in M, 
{ v = f, on (0,T) xdM. 

v G L 2 (0, T; H 2 (M)) D # 2 (0, T; L 2 (.M)) . 



MIl 2 (0,T;H 2 (A4)) ^ ^lu II H; 



(7.3) 



(7.4) 



Estimate (17.41 ) is essentially known, but we give the proof for the readers' conve- 
nience. Let v\ = d 2 v. Then, by hyperbolic estimates, we have 



MIl 2 (0,T;L 2 (.M)) 



<c 



Hi ■ 



On the other hand, since A g v = v±, by the elliptic regularity, we get 



\L 2 (0,T;H 2 (M)) 



<c(lH 



L 2 (Q,T;L 2 (M)) 



+ 



Hi 



Thus, we get 



dv 



di 



<C 



Hz 



Hi • 



Now, for q G L°°(7W), let u> solve 

f {d?-A g +q(x))w 



-q(x)v, in (0,T) x M, 
w{0,-) = 0, d t w(0,-) =0 in M, 

w = 0, on(0,T)xdM, 



(7.5) 
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we apply Lemma [J~T1 to dtw, we can prove 

W W \\l 2 (0,T;H 2 (M)) - C \\qv\\ H l(p tT . L 2( M fi < C\\j\\ Hl 

Thus, for u = v + w, we have 

' (d 2 - A g + q(x)) u = 0, in (0,T)xM, 

«(0,-) =0, d t u(0,-) =0 in At, 

{ u = f, on (0,T) x &M, 

where 

u G L 2 (0,T;H 2 (M))nH 2 (0,T;L 2 (M)), 

and 



(7.6) 



l|A{, g (/)lk 2 <C||/|| Wl . (7.7) 

We shall denote by ||Ag j g||£( Wl ^ 2 ) the operator norm of A gi9 . 

Lemma 7.4 Let (Ai ,g) be a simple Riemannian compact manifold with boundary 
of dimension n > 2, let T > Diam g (.M), there exist constants C > and k G 
(0, 1) swc/j that for any real valued potentials gi, 52 £ =S(Mo) smc/2 f/iaf gi = 92 
on ?«e boundary dAA, we have 

\\Qi-Q2\\ LHM) <C\\Al qi -Al q2 \r c(HuH2) (7.8) 

where C depends on Ai, T, Mq, n, and s. 

Proof . As in (15.111 ). we have 

- H A g,9i ~ A g,g 2 I^CHi,W2) WhWui II5 , a|Il2 ((0)T)x6IA4) 



JdM 



<cx 



2n+5 



l a l||« ll a 2| 



** ll^g, 01 ^s 



^g, q2 \\C(H 1 ,H 2 )- 



(7-9) 



Where 



\H 2n + A {0,T;H 2 (M)) ■ 



Thus, we can complete the proof of (17.81 ) in the same way as in section 5.2. □ 

We set P(z) = (n gi(?1 (z) — n g:92 («)), from Taylor's formula, we deduce for 
1 < J ' < n, and z < 



pW( )=x; 



ip-j 



p=j 



(p- J)! 



■0 (_ T \n-j 

P (p) {z)+ / i — ^-P( n+1 )(r)dr. 



(7.10) 
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Lemma 7.5 There exist C > and jX\ G (0, 1) such that the following estimate 



p(n+l) ( . 



< Ce^ 1 



(7.11) 



holds true for any z < 0. Here C is a positive constant depending on Mq, Ai and 
|| • || denotes the norm in £(H 3 / 2 (dM); H s {dM)). 



Proof . Let h € H 3 / 2 (dM). It follows from Lemma 177X1 

1 

k> 



p( n + 1 )(z)h = -(n + l)\^2 



\ ( x k, qi ~ z 



in+2 



(h,d v 4> kqi )d v (f> k 



qi/^i/Yk,qi 



1 



+ {n + 1)! J2 7\ 2 ^+2 (h,d u ^ kiq2 )d v ^ k:q2 . 



k> 



\ (Xk,q 2 ~ Z) r 



We split P( n+1 )(z)/i into three terms P { : n+1 \z)h = X x {z)h + l 2 (z)h + I 3 (z)h, 
where 



uz) h = - ( „ + 1)! y: l * +2 - j^ xV+2 



l 2 (z)h = -(n + l)\J2 



/ ^ i l(\ ktqi - Z )n+2 (\ km -z) 
1 



{h,d u 4> kAl )dv<f) k , 



m 



fe>l 



(A M2 -^)-+ 2 



X 3 (^ = -(n + l)!^ 



fe ( A ^ 2 " *) 



(h,d u 4> k , qi - d u cf)k,q 2 }dv(j)k, qi 
_ ,_ 2 (h,d v <j) km )[d v <j) kAl - d v (p km ] 



For I\(z)h, we have 



l^i(^)^ll//i/2 ( aA^) < ( re + !) ! IWIz,2 (aM) 
1 1 



£ 

fc>i 



(A*,,, - z)™+ 2 (A fe)92 - z)™+ 2 



|<9^0fc. 



?2ll_ffi/2(ax) 



• (7.12) 



On the other hand, noting that z < 0, X kq . > 0, j = 1, 2, we see that 



1 1 



(A Ml - z)"+ 2 (A Ma - z)«+ 



< Cmax( t^+3,t^+3 j |Afc, gi - A fe , g2 | 

A fc, 9 i A A,g 2 



< 



c 



A: « 



(n+3) |Afc,gi Afc j(j2 | , 



where we have used estimate (1 1 - 15b - On the other hand, since (see (11.141 ) and 

CU) 

\\9u<f>k,q2 \\H 1 / 2 (dM) — Ck™ , 
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we obtain 



i(z)h\\ H i/2( dM ) 



\li(z 



<C 



/ J 2(n+l) 



Afc 01 — AhnA < Ce 



L 2 (dM) Z^ 2(n+i) l A fc,gi ^fc,g 2 l 
fc>l & 



L 2 (dM) 



(7.13) 



For l2{z)h, we have 



l^2(^)^||fri/a(a>i) 



<C 



L 2 {dM) 



£ 



\fe,gi 



19, 



fc>1 l A fc,92 2) 



(7.14) 



Then Lemma 1731 yields 
Afc 



fc>l ^'^ Zj 

Afe, 9l (Afc i9l + Afc : g 2 ) 



cE 



A:> 



\ (Afe,g 2 - z) 



n+2 



Therefore, we find 



For Z%{z)h, we have 



l 2: 2(^)/l||/fl/2( 9 _ M ) < C 11^111,2(9^) ■ 



<Ce. (7.15) 



(7.16) 



l x 3(2)/i|| H i/2 (9 _ M) < C \\h\\ L 2 {dM) J2 T^+T IW*,«i - d u (t)k, q2 \\ H i 



fc>l ^fe,(J2 



\iif i ) > ■ ■■/* n_n i / 2 (9A / l) 

A j. 



^ C ll^lll,2(aM) 2^ 2(n+i) ll^^fc.32 ~ ^0fc,5i[lffV2(aM) 



fe>l 



< C IHIl, 2 ^) 2^ ~2£ \\d u (pk,q 2 - 9 v (f>k, qi \\ H l 
K n 



/ 2 (dM) 



k>\ 



Therefore 

\\ x z{ z ) h \\mi 2 (dM) ^ Ce \\ h \\^{dM) • ( 7 - 17 ) 

The conclusion follows then from a combination of (17.171 ). (17.161) and (17. 13b - □ 

Proof of Theorem |3j From (17.101) and Lemma 17721 we obtain 



P ij) (0) 



<C[\z\~ 3 ~ +\z 



n-j+l pi 



40 



and then 



P ij) (0) 



< C ( izp^ + ur+V 1 



if \z\ > 1. 



In particular 



P (i) (0) < Cmin (p' 2 ^ + p n+l e^) = Ce^ (7.18) 

s p>l \ J 



where p,^ £ (0, 1). Let M g>q be defined as in Lemma[7_3] We can proceed as in the 
proof of Lemma 1731 to prove 



From identity (17- lb . estimates ( 17.191 ) and (17. 18b . we deduce 



(7.19) 



A 8 -A 8 

g,9l g,<22 



c(Hi,n 2 ) 



< Ce^ 4 , 



provided that e is sufficiently small. To finish, we only need to remark that the 
traces of the geometrical optics solutions constructed in section |4]in fact satisfy 

u l\dMi u 2\dM G %1 
so that in the proof of Theorem [Q the right-hand side of (11.12b may be replaced by 



A» -A» 

g,<?l g,92 



This completes the proof of Theorem [3] 



C(Hi,H 2 ) 



□ 
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